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We propose a method to realize a tp Josephson junction by combining alternating and n parts 
(sub junctions) with an intrinsically non-sinusoidal current-phase relation (CPR). Conditions for 
the realization of the tp ground state are analyzed. It is shown that taking into account the non- 
sinusoidal CPR for a "clean" junction with a ferromagnetic (F) barrier, one can significantly enlarge 
the domain (regime of suitable F-layer thicknesses) of the tp ground state and make the practical 
realization of (p Josephson junctions feasible. Such junctions may also have two different stable 
solutions, such as and tt, and cp, or ip and n. 

PACS numbers: 74.45.+C, 74.50.+r, 74.78.Fk, 75.30.+J 



I. INTRODUCTION 

The interest on Josephson junction (J J) devices with a 
ferromagnetic (F) barrier has been continuously increas- 
ing during the last years 1 ^. The critical current den- 
sity j c of a SFS (S - superconductor) junction exhibits 
damped oscillations as a function of the F-layer thick- 
ness d,F so that the Josephson phase <fi can be or tt in 
the ground stated, tt junctions can be used as (non- 
dischargeable) on-chip 7r-phase batteries for self-biasing 
various electronic circuits in the classical and quantum 
domains, e.g. self-biased RSFQ logio 3 - or flux qubits 4 -^. 
In addition, for quantum circuits self-biasing also decou- 
ples the circuit from the environment and improves deco- 
herence figures, e.g., as in the quiet qubit4&£ In classical 
circuits a phase battery may also substitute the conven- 
tional inductance and substantially reduce the size of an 
elementary cell 8 -. Some of these proposals were already 
realized practically^. 

In this context it is even more interesting to create 
a Josephson junction with an arbitrary phase difference 
p (0 < \tp\ < tt) in the ground state — a so-called tp 
junction^. In addition to providing an arbitrary phase 
bias, such ip junctions have rather interesting physical 
properties such as two critical currents, non-Fraunhofcr 
Josephson current dependence on an external magnetic 
field, half-integer Shapiro steps and an unusual behav- 
ior when embedded in a SQUID loop. In long tp junc- 
tions two types of mobile Josephson vortices carrying 
fractional magnetic flux 4>i < $o and $2 = $0 — $1 
($0 ~ 2.07 x 10 _15 Wb is the magnetic flux quantum) 
may exist, resulting in half integer zero field steps, two 
critical values of magnetic field penetration and other 
unusual properties jii 

To obtain a tp junction the Josephson energy density 







(1) 



To achieve this, the current-phase relation (CPR) j(<f>) 
should be different from the usual sinusoidal one 



j{4>) = h sin</>, 



(2) 



where j is the Josephson (super)current density and j-y is 
the critical current density. It has been shown that a tp 
junction can be realized with a second harmonic in the 
CPR, i.e., 



3 (0) = 3i sm 4> + 32 sin 20, 



if j 2 < andii 



2j2 

ji 



> 1. 
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(4) 



should have a local minimum at <p = tp (0 < \tp\ < tt). 



Note that in Eq. ^ j\ is not anymore the critical current 
density, but just the amplitude of the first harmonic and 
j2 is the amplitude of the second harmonic. The critical 
current density j c is then determined by local maxima of 
j(4>) as shown below. 

A non-sinusoidal CPR is not so exotic for Josephson 
junctions with a ferromagnetic barrier—, as demonstrated 
recently in several experiment a 12 ! 13 . For example, in the 
simplest case of SFS junctions consisting of pure S and 
F metals, the CPR is strongly non-sinusoidal at a tem- 
perature T <C T c , where T c is the critical temperature of 
the S metal. Different scattering mechanisms in the F- 
layer also influence the CPR: usual (non-pair breaking) 
scattering ("dirty" limit) 14 ' 15 ; spin- flip scattering^ and 
scattering of electrons from the s to the d-band 5 -^. In ad- 
dition, the transparency of the interfaces or the presence 
of an extra insulating (I) dielectric layer, like in SIFS or 
SIFIS junctions, influences the CPR as wel l 18 ' 19 . How- 
ever, the theoretical models 2 - ^ 15 ' 16 ' 17 ' 18 that take into 
account different scattering mechanisms give similar re- 
sults: The amplitude of the second harmonic oscillates 
and decreases with dp twice faster than the first one. To 
satisfy (|4]) one is tempted to choose <1f in the vicinity of 
0-7T transition, where | j\ \ <C 1 32 \ ■ Unfortunately, there 
j2 > q£ i 15 i 16 i 17 , which excludes a tp ground state. 
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A technique to create a negative second harmonics ar- 
tificially was proposed recentl y 10 ' 20 . By using one of the 
available 0-ir junction technologies 21 ^, one fabricates a 
Josephson junction with alternating short (ji > 0) and 
11 (ji < 0) regions, which all have the simple CPR |2]). 
Here, an effective second harmonic with negative ampli- 
tude is generated for the averaged phase, which is slowly 
varying on the scale of the facet lengt h 10 ' 20 . Let the alter- 
nating and 7r regions have the facet lengths a and b, re- 
spectively, so that a~b. It is assumed that a, b < A}, Aj, 
where Aj = Aj(j°) and Aj = Aj(j^) are the Josephson 
penetration depth in the corresponding parts, with 



Aj(j c ) 



2(i Trd' \ j c 



(5) 



Here n$d' is the effective inductance per square of the 
junction electrodes. 

The effective second harmonic has a maximum am- 
plitude if a = b ~ A,/, and if the values of the critical 
current densities in corresponding regions jf w — jj, see 
Ref. EHO. The effective second harmonic amplitude is 
large enough only if and j\ are very close by absolute 
values, demanding to choose thicknesses dp, a and dF,b 
with very high precision less than 1 A for usual "dirty" 
SFS junctions. This precision is not achievable technolog- 
ically; therefore, a controllable (^-junction is quite hard 
to realize in this way. 

However, ferromagnetic Josephson junctions already 
have some intrinsic second harmonic. Thus, the question 
is what will be the effective CPR in multifacct junctions 
if one takes into account such an intrinsic second har- 
monic. Usually, j2 < in the range of dp where ji is 
rather large. The idea of the method used here is the 
following. We use two thicknesses dF, a and <iF,b where 
j2{dF,a) < and ji(d F , a ) ~ -ji{d F ,b)- By making a 
step-like F-layer changing between dp, a and di?b, we ef- 
fectively cancel the large first harmonic and create an 
effective negative second harmonic which adds up with 
intrinsic second harmonic. 

The paper is organized as follows. In Sec. |TT] we de- 
scribe the model and derive the equations for the effec- 
tive phase using an averaging procedure over the rapid 
oscillations of the phase on the length scales set by a and 
b . Section IIIII contains the discussion and results of cal- 
culations done for a "clean" SFS junction. Section IIVI 
concludes this work. 



II. MODEL 

Mints and coauthor s 20 ' 23 ' 24 considered 1-dimcnsional 
junctions (along x direction) with a critical current den- 
sity ji (x) which is a random or periodic function chang- 
ing on a length-scale a, b <C X a j = Xj{j\ a ,b)- Buzdin 
and Koshele v 10 ' 11 have found an exact solution of this 
problem if j\(x) in ^ alternates between the constant 



values and j\. Both groups assumed that j2 = 
within each region. Instead, here we assume that the 
CPR j(dF(x), (j>) is non-sinusoidal as a function of <f> 
within each region, and alternates between j(dF, a , 4>) and 
j(dp^,(j)) as a function of x. Then the problem of cal- 
culating the phase distribution along such a non-uniform 
structure can be reduced to the well-known problem of a 
non-linear oscillator behavior—. 

The Josephson current in every region, as an odd func- 
tion of the phase, can be expanded in a series of harmon- 
ics 



j(d F (x),4>) = 22j n (d F (x))smncj), 

n=l 



(6) 



and the problem is solved in a general form. Here j n 
denotes the amplitude of the intrinsic n-th harmonic of 
the current density. 

We assume £f *C a,b < Aj, were, £f is the ferromag- 
netic coherence length, i.e. the characteristic length for 
the critical current density nonuniformity at the region 
of the step-like change of junction properties^. Let us 
consider one period of the structure [—a; b) with F-layer 
thickness 



dp(x) 



dF,a, x £ [-a, 0) 
d F .b, x e [0, b) . 



(7) 



We rewrite (|6]) to separate the average and the oscillating 
parts of the Josephson current, i.e., 

oo 

4>) = X] t 1 + Sn{x)] sin ncf), (8) 

n=l 

where j(x,4>) = j{dF(x),4>). The average value of any 
function f{x) is defined as 



(/) = 



a + b 



f{x)dx. 



(9) 



Then the averaged supercurrent of the n-th harmonic is 

aj n (d F . a ) + bj n (d F ,b) /in x 

\3n) = —r , (10) 

a + b 

and the corresponding oscillating part is 

9n(x) = (j n (x)-{j n ))/(j n ), (11) 

so that (g n ) = by definition. Here j n (x,(f>) = 
j n (d F (x),(/>). 

One can represent the Josephson phase (j){x) as a 
sum of a slow component ip{x), changing on a distance 
~ Aj = Xj((ji)), see Eq. ([5|), and a rapid component 
C(x), changing on a distance ~ a, b , i.e., 



4>(x) = i> + c(x). 



(12) 



Here we assume that the junction is short (I <^ A,/) and 
do not write the x dependence for ip. We also assume 
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that the average of fast phase oscillations is vanishing, 
i.e., 



(0 = o, 

and their amplitude is small 

(ICI)«i. 



(13) 



(14) 



The ground state of the Joscphson junction is deter- 
mined by the Ferrel-Prange equation 



2 d 2 (p _ j{x,<p) 
Aj dx^ 



(15) 



Substituting Eqs. ([8]) and (fl2|) into Eq. (fl"5j) and keeping 
the terms up to first order in ((x) we can obtain equations 
for the rapid phase £: 



d 2 C 

A 2 j— - ^2 Pn9n(x) sin nip = 0. 



dx 2 



(16) 



To obtain the equation for the slow phase ip, we average 
Eq. ([Po]) over the length (a + b) <C Aj and get 



d 2 ip 

A 2 — — - };[(3 n sinmp + n(3 n (g n () cosmp] = 0, (17) 



N 



dx 2 



where (3 n = (j n ) The number of harmonics N, 
that is reasonable to take into account within the given 
approximation, follows from the condition N -C 1. 

We have to find the function £(x) from equation (Til)]) , 
calculate average values (g n (), and substitute them into 
(fT7|) . Since for step-like d F (x) the Josephson current 
© is a step-like function of x, Eq. ([16)) has the form 
d 2 (/dx 2 — const on every interval [— a; 0) and [Q;b). Its 
solution is a parabolic segment. The function C(x) must 
be continuous at x = 0, and at the edges, i.e., it must 
satisfy the boundary condition £(— a) = C,(b). Moreover, 
it should satisfy Eq. ([H]) , 

It is convenient to expand the rapid phase as a 

solution of Eq. (flTj)) , into a series 



C(x) = ^ sinn?/;. 



(18) 



n=l 



From this, the average values are calculated as 
2a5j n Sj k 

Pn{9nCk} = . r | ; n, fc=l,2... 

1 0x> I 

where <5j„ = j n (d F . a ) - j n (d F ,b) and 



24A 2 7 ( a + 6)2| (il )| 



(20) 



By definition ||SJ), Aj = Aj |ii a |, where for brevity 
ha = ji(d Fa ), X.j = Xj(j la ). 



The dependence of the Joscphson current on the slow 
phase, that changes on a large distance of the order of 
A j, follows from equation (fT7|) as 



m = \(h)\ 



(3 n sin nip ■ 



^ n(3 n (g n Cfc) cos nip sin fc^i 

n,fe 

(21) ■ 

The CPR (|2T|) contains contributions of two types: in- 
trinsic harmonics, that are given by the first term of (|2ip 
and the effectively generated ones, that are given by the 
second term. Intrinsic contributions are defined by the 
CPR j6]) of the junction regions, while generated contri- 
butions effectively appear as a result of averaging over 
fast oscillations. 

The amplitudes of the effectively generated harmonics, 
which arc proportional to average values (Til?)) , are largest 
by absolute value if a and \Sj n \ reach their maximum. 
This happens if the lengths of a and b facets have the 
largest possible size which still allows averaging, i.e, 



and 



a = b sa A , 



jn(d F ,a) ~ -jn{d,F,b) 



(22) 



(23) 



Condition (|23|) ensures <C |ji a | and consequently 

A ,/ ^> Aj. It was showr>i£, that even if the equality (|2"2"|) 
holds exactly, condition (TTJ]) is satisfied. In this case 
a = l/96|j la | [c.f. m\- 

In all theoretical models developed up to 
no w 2 i 5 i 15 i 16 i 17 i 18 the second harmonic was usually 
considered to be much smaller than the first one (except 
for the points of the 0-7r transition, where j± — ► 0), with 
an even smaller third harmonic. So, the expression ^ 
can be considered as a Tailor expansion in some small 
parameter. Then, keeping terms of the same order of 
this small parameter, we obtain n + k = N + 1 effective 
harmonics in the CPR of the multifacct 0-n junction. 

The number TV of harmonics, that is reasonable to take 
into account, follows from the estimate of the short-range 
phase iVmax \C(x)\ <C 1. As £(x) has a parabolic form, it 
takes its maximum absolute value at the center of every 
interval —a/2 and 6/2. With the estimate 



,(3 n g n (x) sin nip ~ g x 



(19) one obtains 



\C(-a/2)\ 



ab(a + 2b)\Sj 1 \ 
24A> + &)|( J1 >r 



(24) 



(25) 



This is maximized if a = b = Xj. Since max|<5ji| 
2 \ji a \, we estimate max \((x)\ w 1/8. Thus, as a reason- 
able choice one can take N — 3. So, within the above 
approximations, it is reasonable to consider three intrin- 
sic harmonics with amplitudes J1..3 and four generated 
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harmonics, i.e. 

4 

•^(VO = Jn smmp (26) 

n=l 

with 

•A = (ji)+aSj 1 Sj 2 

J 2 = (j 2 ) - a^j'i + 2a8j 1 8js 

h = (j^-SaShSh (27) 

J 4 = -2a5ia - AaSjtSjs . 

Here J ra denotes the total amplitude of the n-th harmonic 
in the CPR for the averaged phase. These expressions 
reduce to the earlier result o 10 ' 11 if one takes into account 
only the intrinsic first harmonic. The expression for every 
harmonic is a simple sum of the corresponding average 
intrinsic harmonic and the generated part. 



A — > 0, equations become linear and their simple ex- 
ponential solutions yield only a sinusoidal CPR. There- 
fore, the most promising strategy for realization of a ip 
junctction is to use a pure SFS junction with a thin Ni 
layer at low temperature (< 0.1T C ). The model describ- 
ing the Joscphson effect in such junctions was established 
long ago 1 ^. It is based on the solution of the Eilcnbcrgcr 
equations. 

What is the measurable critical current density Jc 
(maximum supercurrent) of a Joscphson junction with 
a non-sinusoidal CPR? It is not anymore |ji = |j(7r/2)| 
as follows from ©, but the local maximum of the ex- 
pression © with respect to <j>. Note that © may allow 
several local extrema in the interval [0, 2tt). 

Figure [T]Ja) shows the three first harmonics vs F-layer 
thickness and Fig. QJb) shows the corresponding Jc(dp) 
for the model of an SFS junction as described in Ref. [Tj 
(see also Ref. [Tol ). It is interesting to note, that near 



III. RESULTS AND DISCUSSION 

We first address the question which type of ferro- 
magnetic junction with 0-7r facets can satisfy the con- 
ditions for a tp junction in the best way. The critical 
current density j\ of a Joscphson junction with a ferro- 
magnetic barrier having a complex coherence length £p 
(£p = S,^ 1 + ) decays oc cxp(— c£f/£i) and oscillates 
with a period 2 - 2tt^2 in the "dirty" limit (electron mean 
free path I <C £1,2)' The decay length £1 depends on 
I as well as on a pair-breaking scattering length in the 
ferromagnet (spin-flip scattering^ or scattering into the 
d-band 2 ^.; both influence the CPR). Pair-breaking scat- 
tering leads to 17 ' 28 £1 < £2 and the Joscphson current 
decays so rapidly in the "dirty" limit, that there are 
only very tiny regions of dp, a and dpj, on the dependence 
ji{dp) where ji(dp ya ) ~ —ji(dp,b)- In the "clean" limit 
the Josephson current decays much more slowly with in- 
creasing dp. The cleaner is the ferromagnet, the larger is 
I, and the slower is the decay; c.f. Eqs. (23) and (24) with 
Eq. (7) from Ref. [H, or see the discussion in Ref. lU In 
the limit I 3> (dp 1 , £1,2) the critical current density de- 
creases as 1 ^ 1/dp. Moreover, a "clean" SFS junction has 
a non-sinusoidal CPR, and its second harmonics j-2 < 
in some regions of dp far from 0-7r transitions, that can 
effectively help to satisfy the conditions for the realiza- 
tion of a ip junction. It was shown in different models, 
that the second harmonic ^(^f) decays and oscillates 
with dp twice faster than j\(dp). Therefore, it is reason- 
able to take an F-layer with dp.a and dp^ of the order 
of a few £2- Usually, £2 = vf/2E cx , where vp is the 
Fermi velocity and E cx is the exchange magnetic energy 
in the ferromagne t 15 ' 27 i 29 . It was established experimen- 
tally, that for pure ferromagnetic metals £2 is largest for 
N i 30 ' 31 ' 32 i 33 . Higher harmonics also decay rapidly with 
increasing temperatur o 15 i 34 . We have also checked this 
statement for models described i n 14 i 17 . It is clear that 
by approaching T c from below the superconducting gap 
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FIG. 1: (Color online) Properties of a uniform SFS Joseph- 
son junction in the "clean" limit at T = 0.1T C : (a) ampli- 
tudes of the first three harmonics ji, ]2, J3 of Eq. ([6j| as a 
function of the ferromagnet thickness d_F- (b) measurable 
critical current densities Jci and \Jm\ vs dj?, which are re- 
alized at the phase < 4> < n and tt < (j> < 27r, respec- 
tively. Inset shows for dF = 0.5^2 the CPR j(cj>) from (3) 
in Ref. [l^ (dashed line) and from the approximation by the 
first 3 harmonics (dotted line), with the corresponding Ej((j>) 
dependence (solid line). For comparison, in (b) the critical 
current |ji(<Zf)| of a junction with only the first harmonic in 
the CPR is also shown. Here, pairs of points (connected by 
arcs) where ji(dp t a) = —ji(dF.b) correspond to a (p ground 
state of a multifacet SFS junction with alternating thicknesses 
dF, a , dF,b (see areas 1,2,3 in Fig. HJa)). 
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a 0-7r transition, when the first harmonic is small, j(tp) 
may have two different local maxima, and the measurable 
Jc(dp) has two different values Jci, \Jc2\ depending on 
the initial state of the junction^ (see Fig. Hfb)). An 
example of a CPR j(4>) with two different maxima and 
the corresponding Josephson energy Ej ((f)) ^ near a 0-tt 
transition are presented in Fig. [IJb) (inset). Thus, even 
for a uniform Josephson junction (where the (^-ground 
state is impossible) close to a 0-tt transition, two different 
values of critical current density could be realized if its 
CPR differs enough from the sinusoidal one. 

In the framework of the clean SFS junction model^ 
we investigate a multifacet junction with a = b = Xj 
A,/. For A,/ we take the value corresponding to the first 
minimum of ji at <1f ~ £2, see Fig. HJa). For clean SFS 
junctions this value is 3 - 2 - w 10 kA/cm 2 , which corresponds 
to Xj w 3 fim and which allows to realize a w b with 
reasonable precision. Below we investigate the ground 
states in a Josephson junction with alternating regions of 
length a, b and F-layer thicknesses d,F, a an d dF,b varying 
from to a few £ 2 . 

The ground state corresponds to a local minimum of 
the energy (JTJ with the CPR (|2ip. The junction has a 
stable static solution ip = (0-phasc) if 



has at least one real solution z = cos tp, satisfying the 
conditions |z| < 1 (which gives < \ip\ < tt), and 



N 

E 

n=l 



W+l-ri 



0") ~ a E nk5 3n5jk 



fc=l 



> 0, 



and the solution ip = tt (?r-phase) if 



N 

E 

n=l 



JV+l-n 



(i„)-a (-l) n+k nkdj n 5j k 



> 



(28) 



(29) 



For the CPR ([26]) these conditions have the following 
form: For the 0-phasc 

Ji + 2 J 2 + 3 J 3 + 4J 4 > 0, (30) 

and for the 7r-phasc 

Ji - 2J 2 + 3 J 3 - 4J 4 < 0. (31) 

Both solutions -0 = and ip = tt coexist when these 
conditions are satisfied simultaneously, i.e. 

2J 2 +4J 4 > \Ji+3J 3 \. (32) 

If both conditions (|30"|) and (|3Tj) are not satisfied, i.e. 

2J 2 + 4J 4 < -|Ji+3J 3 |, (33) 

only the 99 ground state is possible. The conditions (|30p 
(|3"3"|) coincide with the conditions Q if J3 = J4 = 0. 
Generally, the tp junction is realized if the Josephson en- 
ergy Ej(ip) ((T|) has a local minimum, i.e., if the CPR 
J(ip) (|26|) crosses J = from a negative to a positive 
value at some point ip = tp ^ 0, n. This is the case if the 
equation 

J(ip) - 8 J 4 z 3 + 4 J 3 z 2 + (2J 2 - 4J 4 )z + Ji- J 3 = 0, (34) 



5J 
9^ 



~ 32 J 4 z 4 + 12 J 3 z 3 + (4 J 2 - 8 J 4 )z 2 + 
+ (J x -9J 3 )z-2J 2 + 4J 4 > 0, 



(35) 



that ensures the local minimum of Ej(ip). Starting 
from a pair of dp^ a , d,F,b we calculate J1...4. Then from 
Eqs. (j30|) - ([35)) the possible ground states are identified. 
The resulting phase diagram is shown in Fig. [21 where dif- 
ferent ground states for each pair of (If at d,F,b are shown 
by different colors. 

Figure [2ja) shows the results obtained if only the first 
intrinsic harmonic is taken into account. Here, the areas 
of and tt ground state phase are separated by slim re- 
gions of tp phase. It is clear that the phase diagram is 
symmetric with respect to the line dp, a = dp,b- There- 
fore, below, without loosing generality, we focus on the 
case dF.a > dp.t- In the chosen interval of thicknesses 
there are 3 areas of tp phase, marked as 1,2 and 3. In 
these areas ji(dF,a) ~ —ji(dF,b), which corresponds to 
the pairs of dp,a and G?i? & shown in Fig. [ljb) as 1,2 and 
3. If dp^ a and dp^ are not well controllable, the area S^p 
of tp regions is proportional to the probability of the tp 
junction realization. For all three different cases shown in 
Fig-[2Ia) S v is rather small (< 0.025£|). Hence, one has 
to control cZf.o and dF,b extremely precisely in order to re- 
alize a tp junction. In the case of a "dirty" SFS junction, 
when j'i(cZf) decays exponentially, the area S v is even 
smaller, and the probability to fabricate a (^-junction is 
vanishing. 

However, if we take into account a non-sinusoidal CPR 
for the "clean" SFS junction, the areas of the tp ground 
state become much larger, as can be seen in Fig. HJb). 
This is a consequence of the fact that the intrinsic sec- 
ond harmonic j 2 (c?_f) < in the corresponding regions, 
which efficiently helps to make the absolute value of the 
generated second harmonic large enough. 

In particular, the areas 2 and 3 from Fig. Ufa) merge 
and form a compact tp ground state region with area 
S v w 0.6£| around dF, a ~ 2.8£ 2 and d F ,b ~ l-4£ 2 . This 
region seems to be very well suited for the experimental 
realization of a tp junction, as it does not demand to pro- 
duce an extremely thin F-layer and at the same time al- 
lows for reasonably large tolerances in sample fabrication. 
Taking the ferromagnetic coherence length £ 2 ~ 1.2 nm 
from experiment a 31 ' 32 , the linear size of the tp region in 
FigEJb) is Ad F , a ~ Ad Fi6 - 0.5 ... 1 nm. Modern tech- 
nology allows the control of dp with such precisio n 35 ! 36 ' 37 . 

Finally we note that the structure considered here may 
have simultaneously two different stable static solutions: 
and tt, and tp, or tt and tp, as can be seen in Fig.^b). 
The corresponding Ej(ip) curves are presented in Fig. [3l 
where one can see several local minima of the Josephson 
energy as a function of the phase. Probably, taking into 
account the next harmonics, it is also possible to find 
two different tp ground states. However, we expect the 
corresponding regions within the dp ta ,dF.b plane to be 







1 2 3 

ferromagnet thickness d /1L 



12 3 

ferromagnet thickness d p it 



FIG. 2: (Color online) Ground states scheme for the multifacet SFS junction with F-layer thickness periodically changing 
between di?a and dp,b plotted on the (dF, a , dF,b) plane. The values of the ground state phase corresponding to each region 
are shown, and the areas for some of the different ip ground state regions are also indicated. In (a) the calculation includes 
only the first harmonics in every uniform part, i.e., two generated harmonics. The ip- regions 1,2 and 3 in (a) correspond to 
the pairs of dF,a and dF,b that are shown by arcs in Fig. QJb). In (b) the first 3 harmonics and the corresponding 4 generated 
harmonics are taken into account, see Eqs. (|26[) and (|27[) . The dependence Ej(ip) for points A,B,C,D is shown in Fig. [3] 
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FIG. 3: (Color online) Josephson energy Ej(ip) of multifacet 
SFS junctions for different combinations of F-layer thicknesses 
dF, a and dF,b as indicated by points A,B,C, and D in Fig. Ob). 
A and D correspond to the tp ground state; B corresponds to 
two stable states ir and ip; and C corresponds to two stable 
states and ip. 



so small that it will be quite hard to fabricate such a 
structure. 



IV. CONCLUSION 



We propose a method to realize a ip junction, i.e. with 
the Josephson phase (p ^ or 7T in the ground state, 
based on SFS junctions. In a uniform SFS Josephson 
junction near a 0-7r transition (with appropriate F-layer 
thickness (If) the first harmonic ji — > 0, while the sec- 
ond harmonic j% dominates, however usually with posi- 
tive sign (j*2 > 0), which excludes the formation of a ip 
ground state. Instead, the main idea of the method used 
here to make a ip junction is the following: We choose a 
thickness cLf where the second harmonic is large and neg- 
ative. To cancel the first harmonic we use a periodic step- 
like modulation of the F-layer thickness between <ii? a and 
d,F,b- Here, d,F, a and d,F,b are chosen such that for both 
of them the second harmonic is negative, while the first 
harmonic is j\(dF,a) ~ ~ .7i(^f,&)- Periodic modulation 
not only cancels the first harmonic but also generates 
an additional negative second harmonic for the average 
phase. This effect can be made stronger if one works with 
a "clean" ferromagnetic barrier and uses alternating re- 
gions of equal length a w b ~ Aj. 

Different mechanisms, leading to a significantly non- 
sinusoidal current-phase relation of SFS junctions are an- 
alyzed. The CPR mostly different from the sinusoidal one 
is obtained for a "clean" SFS junction at low tempera- 
ture. In this case, there are reasonably large regions of 
thicknesses d,F, a and dF,b (m comparison with the descrip- 
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tion taking into account only the first harmonic) where 
the multifacet Josephson junction has a (^-ground state. 
Moreover, for some values of dp, a and d,F,b such structures 
may have two different ground states (two local minima 
of the Josephson energy as a function of the phase) : and 
7r, and ip, or 7r and ip. Our analysis gives some practi- 
cal recommendations for the fabrication of SFS junctions 
with arbitrary phase shifts tp in the ground state. 
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